Renormalized One-loop Theory of Correlations in Disordered Diblock Copolymers 
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A renormalized one-loop theory (ROL) ' is used to calculate coiTections to the random phase approximation 
(RPA) for the structure factor S{q) in disordered diblock copolymer melts. Predictions are given for the peak 
intensity S{q*), peak position q*, and single-chain statistics for symmetric and asymmetric copolymers as 
functions of where x is the Flory-Huggins interaction parameter and A'' is the degree of polymerization. 
The ROL and Fredrickson-Helfand (FH) theories are found to yield asymptotically equivalent results for the 
dependence of the peak intensity S{q*) upon x^ for symmetric diblock copolymers in the limit of strong 
scattering, or large x-^> but yield qualitatively different predictions for symmetric copolymers far from the 
ODT and for asymmetric copolymers. The ROL theory predicts a suppression of S{q*) and a decrease of q* 
for large values of xN, relative to the RPA predictions, but an enhancement of S{q*) and an increase in q* 
for small X-^ (X^ < 5). By separating intra- and inter-molecular contributions to S~^{q), we show that the 
decrease in q* near the ODT is caused by the q-dependence of the intermolecular direct correlation function, 
and is unrelated to any change in single-chain statistics, but that the increase in q* at small values of xN is a 
result of non-Gaussian single-chain statistics. 
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I. INTRODUCTION 

Disordered diblock copolymer melts exhibit composition 
fluctuations that can be measured by small angle X-ray 
(SAXS) and neutron (SANS) scattering experiments. Results 
of such experiments are often analyzed by fitting the scat- 
tering intensity as a function of wavenumber q to Leibler's 
random-phase approximation (RPA) theory for the structure 
factor S{q)? The RPA provides a rather accurate description 
of S{q) in melts of very long copolymers far from the order- 
disorder transition (ODT), but fails near the ODT of nearly 
symmetric diblock copolymer, where strong composition fluc- 
tuations cause a breakdown of the underlying self-consistent 
field (SCF) approximationji 

The scattering intensity measured by small angle scatter- 
ing from AB diblock copolymer melts is proportional to the 
structure function S{q) = / (ir((5c^(r)(5cyi(0))e''i '", where 
^C/i(r) represents a deviation of the number concentration of 
A monomers from its spatial average, and q = |q|. Scattering 
from a disordered diblock copolymer melt generally exhibits 
a maximum S{q*) at a nonzero wavenumber q*. 

Leibler's RPA theory predicts an inverse structure function 
of the form 



cNSo\q) = Fo{qR)-2xN. 



(1) 



Here, x is an effective interaction parameter, N is the degree 
of polymerization, and c is the monomer number concentra- 
tion. The dimensionless function Fq {qR) has a minimum at a 
wavenumber qo = xq /R, yielding a corresponding maximum 
in 50(9), where R cx ^/N is the copolymer radius of gyra- 
tion and is a dimensionless number that depends upon the 
copolymer composition, but that does not depend on x- Here 
and hereafter, we use 5o((z) to denote the RPA approximation 
for S{ci), and go to denote the RPA prediction for q* . 

Eq. ([T]l has been found to adequately describe the q- 
dependence of most SANS and SAXS experiments,"^-^ and has 
been widely used to extract values for the interaction param- 
eter x(r) as a function of temperature T.^'^ Careful compar- 



isons of eq. ([U to both experimental and simulation results 
have, however, also revealed some limitations, particularly 
near the ODT: 

The RPA predicts an inverse peak intensity 



cNS^\qo)^2[{xN)s-xN] 



(2) 



that depends linearly on x- Here {x^)s — FQ{qQR)/2 is 
the mean-field spinodal value. For symmetric diblock copoly- 
mers, (x-^)s = 10.495. To make a meaningful comparison of 
the RPA to a SANS or SAXS experiment, x must be allowed 
to be a function of temperature T, which is usually fit to the 
form x(r) ~ A/T + B, but is assumed to be independent 
of chain length N and wavenumber q. Sufficiently far away 
from the ODT, the dependence of SANS and SAXS data on T 
and N is often adequately described by the resulting theory. 
In data taken near the ODT, however, plots of the inverse peak 
scattering intensity vs. 1/T exhibit a characteristic nonlinear- 
ity that cannot be described by the RPA.^'^""' 

The RPA predicts a peak wavenumber q* = qo that is 
proportional to the inverse radius of gyration of a Gaussian 
chain. The peak wavenumber q* is thus expected to vary with 
chain length iV as cx N^^^"^ at fixed temperature. Stud- 
ies of homologous series of diblocks of differing N at fixed 
T have found q cx N~^/^ far from the ODT, but a signifi- 
cantly stronger A^-dependence near the ODT.'^ '** In any sin- 
gle sample, the RPA suggests that qo should depend on tem- 
perature only as a result of the slight intrinsic temperature de- 
pendence of the monomer statistical segment lengths. Numer- 
ous experiments have shown that q* decreases with decreas- 
ing temperature at a rate that systematically increases near the 
0^-^^5,6^-23^5,16^-21. jjjj^j ^Yicit appears to be too large to be 
explained by the observed temperature dependence of the pure 
component statistical segments lengths."-'^ 

Simulation studies have provided analogous results for 
both the peak intensity S{q*) and the peak wavenumber q*, 
though for chains shorter than those studied in most exper- 
iments. Both lattice Monte Carlo and molecular dynamics 
simulations have shown a strongly non-linear dependence of 



S~^{q*) on inverse temperature^"— a gradual decrease in 
q* with decreasing temperature,—"— and non-Gaussian chain 
statistics. Though the shift in q* near the ODT was ini- 
tially described as a result of "chain stretching",'^ Binder and 
Fried^** found that the temperature dependence of q* in lat- 
tice Monte Carlo simulations was significantly stronger than 
the temperature dependence of the radius of gyration, and so 
suggested that the shift in q* might be caused primarily by 
changes in intermolecular correlations. Similar results were 
later obtained in an analogous experimental comparison by 
Bartels and Mortensen,^' in which S{q) and Rg were mea- 
sured independently. 

Leibler originally presented the RPA for S{q) in diblock 
copolymer melts as part of the more general analysis of 
the weak-segregation limit of the self-consistent field theory 
(SCFT) for such systems.^ Leibler found that SCFT predicts 
a second-order (continuous) transition between the disordered 
and lamellar phase in the special case of symmetric diblock 
copolymers, and a weakly first order transition for nearly sym- 
metric copolymers. Leibler also noted, however, that Bra- 
zovskii^" had previously analyzed a phenomenological model 
of transitions between a homogeneous liquid and a periodic 
state in systems for which Landau theory predicts a second- 
order transition, and had concluded that such systems actu- 
ally always exhibit a fluctuation-induced first-order transition. 
Brazovskii also presented a theory for the dominant correc- 
tions to the Landau (or RPA) theory for S{q) near such tran- 
sitions, which were found to be unusually large for this class 
of system. It was thus clear from the outset of interest in this 
subject that a more sophisticated treatment of fluctuation ef- 
fects would be needed to adequately describe the vicinity of 
the ODT. 



A. Fredrickson-Helfand Theory 

The task of incorporating fluctuation effects into Leibler' s 
theory was undertaken by Fredrickson and Helfand.^ 
Fredrickson and Helfand approached the problem by con- 
structing an approximate mapping of Leibler's theory for 
block copolymer melts onto Brazovskii's theory of weakly 
first order crystallization. The Brazovskii and FH theories 
are both based upon an expression for the partition function Z 
as a functional integral of a composition order parameter field 
ipir), of the form 

Z = J D[i:] e-^W/*^--^. (3) 

Fredrickson and Helfand began their analysis by assuming 
that the effective Hamiltonian in eq. (O can be approximated 
by the SCFT free energy functional, while using Leibler's 
Taylor expansion of this quantity. We will refer to this as- 
sumption as a mean-field effective Hamiltonian (MFEH) ap- 
proximation. They then introduced a variety of further mathe- 
matical approximations in order to map the SCFT free energy 
functional to the relatively simple phenomenological expres- 
sion considered by Brazovskii. 



The FH theory yields a prediction for the inverse structure 
factor S^^{q) as a sum 

S-\q) = S^\q)+SS~\q), (4) 

where So{q) is the RPA structure factor, eq. ([T]i, and where 
6S^^{q) is given by a self-consistent equation 

cN5S-\q) = J^^, (5) 

in which 

T = cNS-\q*) (6) 
is a normalized inverse peak intensity, and 

N = N{cb^f (7) 

is the so-called invariant degree of polymerization. 

The parameter N^^^ that appears in eq. (|5]l is a measure of 
chain overlap: 7V^/^ is proportional to the number of chains, 
each occupying a volume N/c, that can fit in the volume 
i?3 ^ N^/^b^ explored by any one chain. Values of N in ex- 
perimental studies on diblock copolymer melts near the ODT 
have rang ed from about 500 (ref.'^) to 10000 (ref."), corre- 
sponding to systems with molecular weights that range ap- 
proximately from 10 to 200 kg/mol. 

The mathematical approximations used by Fredrickson and 
Helfand yield an expression for SS^^{q) that is independent 
of q, as indicated in eq. (|5]l. This yields a peak wavenumber 
q* equal to that predicted by the RPA. A subsequent extension 
of the theory by Barrat and Fredrickson^' (BF) allowed for the 
possibility of a fluctuation-induced shift in q*, and predicted a 
peak wavenumber q* that decreases with increasing xiV near 
the ODT, as seen in experiments and simulations. 

Several experimental studies have quantitatively compared 
data for both the peak intensity and S{q*) and the peak 
wavenumber q* for model diblock copolymer melts near the 
ODT to the FH and BF theories. The FH theory has been 
found to describe the temperature dependence of the peak in- 
tensity near the ODT reasonably well in several systemsj^i^ 
Almdal et al}'' have also compared the N dependence of q* 
in a series of symmetric diblock copolymers at constant tem- 
perature to the predictions of the BF theory, and also reported 
reasonable quantitative agreement. 

Notwithstanding its success in describing many aspects of 
the experimental results, the FH theory has several shortcom- 
ings that are inherent in how it was derived^^: 

(1) The use of the SCF free energy functional as an effec- 
tive Hamiltonian (the MFEH approximation) has no rigorous 
basis. 

(2) Predictions of both the full one-loop MFEH theory and 
the approximation studied by Fredrickson and Helfand are 
very sensitive to the effects of short-wavelength (monomer 
scale) fluctuations that these coarse-grained theories cannot 
accurately describe. In the jargon of field theory, the theory is 
ultraviolet (UV) divergent. 

(3) Mathematical simplifications introduced by Fredrickson 
and Helfand limit the potential range of validity of the theory 
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to wavenumbers q q* at temperatures very near the ODT in 
melts of long, nearly symmetric diblock copolymers. 

The UV divergence of the FH theory is a complication that 
was not necessarily fatal, but that was not (we think) taken 
sufficiently seriously in early work on this subject. Fredrick- 
son and Helfand simply ignored all UV divergent contribu- 
tions to their expression to SS~^{q) without commenting on 
the existence or possible physical interpretation of the diver- 
gence, and simply reported the UV convergent parts of the 
resulting integrals. Such a procedure can generally be justi- 
fied if and only if it can be shown that the neglected terms 
can be absorbed into a renormalization of the values of few 
phenomenological parameters, such as the x parameter Kud- 
lay and Stepanow^^ later analyzed the UV divergent contri- 
butions in the MFEH theory for SS^^{q), and asked whether 
they could be absorbed into a simple renormalization of the 
value of the interaction parameter but concluded that this 
interpretation was untenable. 

B. Renormalized Auxiliary Field Theory 

More recently, several authors have contributed to the de- 
velopment of a set of very closely related renormalized aux- 
iliary field theory of corrections to the RP A*'^^-^^ — that has 
both a more rigorous theoretical basis and a potentially wider 
range of validity than the FH and other MFEH theories. We 
will refer to this as a renormalized one-loop (ROL) theory. 
The present paper presents predictions of the ROL theory de- 
veloped in ref. ' for correlations in disordered diblock copoly- 
mer melts. 

The way that the ROL theory is derived avoids most of the 
aforementioned limitations of the FH theory: 

(1) The theory has a rigorous starting point: It is based 
on the Edwards auxiliary field representation of the partition 
function Z for a simple coarse-grained model. 

(2) It has been shown' that the UV divergence of the one- 
loop auxiliary field theory for S{q) can be removed by a renor- 
malization procedure in which all contributions that are sen- 
sitive to monomer scale structures are absorbed into shifts in 
the values of a few phenomenological parameters, i.e., of the 
effective interaction parameter x and the monomer statistical 
segment lengths. 

(3) The range of validity of the theory is not intrinsically 
limited to the vicinity of the ODT, or to wavenumbers near 
q*. 

The ROL theory yields a prediction for cNS ^{q) as the 
sum of an RPA contribution of the form given in eq. ([T), with 
renormalized values for the x and b parameters, plus a univer- 
sal correction of the form 

cNSS-\q) = .^H{qR,xN,fA,bA/bB). (8) 

No simple analytic form exists for the dimensionless function 
H, which is evaluated here by numerically integrating a set of 
related Fourier integrals (see Sec. Ill and App. A for details). 

We have argued previously ''^^ that the ROL theory appears 
to be the first correction to the RPA in a systematic expan- 



sion of an underlying universal expression for cNS ^ (q) in 
diblock copolymer melts as a function 

cNS-\q) = D{qR,xNjA,bA/bB,N) (9) 

that depends only on the parameters that appear in the RPA 
and on N. It appears that the ROL theory is part of an ex- 
pansion of this underlying function in powers of N^^^'^, in 
which an RPA theory with renormalized parameters is recov- 
ered in the limit N oo, and in which the ROL theory is the 
dominant 0{N^^^^) correction. This claim is based upon a 
simple power counting argument^^ that shows that, if all di- 
agrammatic contributions to S^^{q) within an infinite series 
expansion^^ are organized into a loop expansion, the UV con- 
vergent part of the L-loop contribution will be multiplied by 
a prefactor of N^^/"^. This argument assumes, however, that 
a renormalization procedure similar to what has been used to 
interpret and remove all UV divergent contributions from the 
one-loop theory for S{q) can be extended to higher-order con- 
tributions to the loop expansion. If this is so, the accuracy of 
the one-loop theory should increase systematically with in- 
creasing N , for any value of x^, and for asymmetric as well 
as symmetric copolymers. 

We have also discussed the mathematical relationship be- 
tween the ROL and FH predictions for S{q) in more detail 
elsewhere. We showed that, despite differences in how these 
two theories were derived, they yield asymptotically equiv- 
alent results for the divergence of 5S^^{q*) for symmetric 
block copolymers near the ODT, where S{q*) is large. We 
found that the asymptotic behavior of the ROL prediction for 
SS{q*,xN) for symmetric copolymers at g = go near the 
spinodal xA^ is given by eq. (|5]l, with the same value for the 
numerical prefactor B as that found by FH. The FH theory is 
thus a correct asymptotic approximation to the full ROL the- 
ory for symmetric copolymers, valid for long polymers suffi- 
ciently near the ODT. Significant differences between the pre- 
dictions of the two theories are thus possible only for symmet- 
ric diblock further from the ODT, and for asymmetric copoly- 
mers. 



II. CORRELATION FUNCTIONS 

We consider a melt of AB diblock copolymers in which 
each chain contains N monomers, or f iN monomers of type 
i, where i = A or B. Let c = be the overall monomer con- 
centration and p — c/N be the concentration of molecules. 
Let bi and k ~ v/bf denote the statistical segment lengths 
and packing lengths for monomers of type i, respectively. The 
fluctuating local concentration Ci{r) of i monomers is given 
by a sum Ci(r) = J2,n,s^i'^ ^ R„,i(s)), in which R„„;(s) is 
the position of monomer s on block i of molecule to. 

We are primarily interested in the behavior of the structure 
factor matrix (in this section, the theory is presented using 
wavevectors q as arguments) 

5,, (q) ^ fdr (,5c,(r)<5c, (0))e''^•^ (10) 
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where 5ci(r) = Ci(r) — (q) is the deviation of the monomer 
concentration Ci(r) from its mean value (q) = cfi. Anal- 
ogously, we also define an intramolecular correlation func- 
tion, ily (q), that arise from correlations between pairs of 
monomers on the same chain. 

To distinguish the effects of intra- and inter-molecular cor- 
relations, it is useful to introduce a generalized Omstein- 
Zemicke equation '•^''^^ 

Sr^\q)^nT^\q)-a,ici), (11) 

in which ily (q) denotes the true intramolecular correlation 
function (rather than the random walk approximation used by 
Leibler), and (q) is a direct correlation function that is de- 
fined by eq. (fTTT l. 

In a dense, nearly-incompressible liquid with monomers of 
equal volume, we may assume that the eigenmodes of the 2x2 
matrix Sij (q) for values of q of order the inverse coil size are 
given approximately by a fluctuating composition mode, of 
the form ((5c^(q), (5cB(q)) c>c (1,-1), and a density fluctua- 
tion mode, of the form ((5cA(q), <5cB(q)) oc (1, !)• Fluctua- 
tions of the monomer density mode are strongly suppressed by 
the low compressibility of the liquid, while composition fluc- 
tuations can be quite large. In the limit of negligible density 
fluctuations, one can define a scalar structure function 

S{q) = Saa{ci) = SBB{q) = -SAB{q). (12) 

By combining this assumption of incompressibility with the 
Ornstein-Zemicke equation, it is straightforward to showiiii*^ 
that S*^^ (q) may be expressed in the generalized RPA form 

cNS-\q)=F{q)-2xa{q)N, (13) 

in which the functions 

^^(q) = ciV^f]r.i(q)e,e^., (14) 

are defined by projecting il^^(q) and — Cij (q) onto the sub- 
space of pure composition fluctuations, where e = {eaj^b) = 
(1,-1). Written more explicitly, this yields expressions 

F(q) = cNn{q)/W{q), (16) 

Xa(q) = ^ [Caa{ci) + CBB(q) - 2CABiq)] , (17) 

in which 

f^(q)EE^a,(q), 

ij 

W{q) = nAA{<l)nBB{q) ~ ^ab{^)- (18) 

We will refer to the wavevector-dependent function Xa(q) de- 
fined above as the "apparent" interaction parameter 

Despite a superficial similarity, eq. (fTjt is not equivalent 
to the RPA approximation, eq. ([TJ. Eq. (fT3] l simply define 
the functions F{q) and Xa{q), by relating these quantities to 
the correlations functions S{q) and without assuming 



anything about the behavior of these functions. The RPA is 
instead obtained by approximating the intramolecular corre- 
lation functions used to define F{q) by those of a Gaussian 
chain, and approximating Xa (q) by a parameter x that is in- 
dependent of both wavenumber q and chain length N. The 
most general form of eq. ([U can allow x to exhibit an arbi- 
trary dependence upon both temperature and composition, but 
would cease to have any predictive value if it also allowed for 
an arbitrary dependence on q and N. 



III. RENORMALIZED ONE-LOOP THEORY 

In this section, we briefly review the main results of the 
ROL theory for S'(q) in diblock copolymer melts.— 

The derivation of the ROL theory is based on a coarse- 
grained model in which the total potential energy U is the 
sum of intramolecular potential C/chain of a set of Gaussian 
chains plus a short-range non-bonded pair interaction. We 
consider a pair potential of the form Uij{r) = uoeyi^ACr), 
where (r) denotes a short-range function with a unit inte- 
gral, J dr6A{r) = 1, and with a characteristic range of in- 
teraction A^^. The interaction parameters eij are given by 
eAA = £BB = Bo and e^s = Bq + xo, where Bq is a di- 
mensionless compression modulus, and where xo is a "bare" 
Flory-Huggins interaction parameter 

The one-loop approximation for Sij (q) given in ref. ' was 
based on a more general diagrammatic expansion that was 
discussed in refi^. There, the details were shown on the 
construction of expansions of various correlation functions in 
terms of cluster diagrams that represent Gaussian chains in- 
teracting via a screened interaction G. The Fourier transform 
G(k) of this screened interaction is given by 

Gij\k) = n{k) + u-^{k), (19) 

in which (k), t/y (k) and r2ij(k) are 2x2 matrices, and 
inversion refers to matrix inversion. This is a straightforward 
generalization to multi-component systems of the screened 
interaction introduced by Edwards to describe excluded vol- 
ume interactions in concentrated homopolymer solutions. The 
limit of incompressible dense liquids is obtained in this for- 
malism by taking Bq oo. 

The one-loop approximation for S{q) is based upon an 
analysis of the two diagrams shown in Fig. 1. The diagram 
on the left yields the one-loop contribution to the intramolec- 
ular correlation function riy (q), while the one on the right 
yields a one-loop contribution to the direct correlation func- 
tion (q). 



A. Intramolecular Correlations 

The one-loop theory predicts an intramolecular correlation 
function 

n,j{q)=n,j{q)+An,j{q), (20) 
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in which fjy (q) is the correlation function for a gas of non- 
interacting Gaussian chains. The correction Aily (q) arising 
from interactions is given by a Fourier integral 

Af],,(q) = -^^7/^g,(q,-q,k,-k)G«(k). (21) 

Here and hereafter, summation over repeated subscripts is 
implicit, and = {2tt)~^ J dk. The four-point function 

Tpljli ('^' is defined by 

^^jki = figL(q, -q, k, -k) - %(q)6fe,(k), (22) 

in which fiA:i(k) = flki(k)/p is a single-chain correlation 
function, normalized by the molecular density p — c/N, and 
fi^^],;(q, — q, k, — k) is an analogous single-chain four point 

function for monomers within blocks i, j, k, and I. ' The phys- 
ical content of eq. (ISTT i is shown schematically in the left di- 
agram of Fig. [T] in which the curve represents a single chain, 
along which two segments interact via a screened potential G. 




qji jf-q 



FIG. 1 . Diagrams illustrating the effects of the screened intra- (left) 
and inter- (right) molecular interactions on the correlations of seg- 
ments of type i and j. Additional segmental indices fc, I, m, n runs 
through the whole chain, k and q are wave vectors for the interac- 
tions or correlations, and k± = k ± ^ is one particular choice made 
to conserve momentum. 



The integral in eq. (ISTT i is ultraviolet (UV) divergent: The 
value of the Fourier integral in this equation is dominated by 
contributions of large wavevectors, and diverges in the ab- 
sence of a high-k cutoff. If the range of k is restricted to 
|k| < A, the value of the integral increases linearly with the 
cutoff wavenumber A for KRg ^ 1. 

It was shown in the previous work,' however, that this 
type of sensitivity of the one-loop prediction for ^ij(k) to 
short-wavelength correlations can be entirely attributed to the 
change in the value of the effective statistical segment length. 
Specifically, it was shown that the ROL result for fii^ (q; A) 
with a cutoff wavenumber A is consistent with an expression 
of the form 



n,j{q) = n,j{k; b{A)) + Sn,j{q), 



(23) 



in which ilij (k; 6) is the correlation function for a Gaussian 
chain with a cutoff-dependent statistical segment length b{A) 
that is different than that of a polymer in vacuum, as the re- 
sult of interactions between polymers in a dense liquid, and 
in which Sftij (q) is a small cutoff-independent correction to 



Gaussian statistics. The treatment of intramolecular correla- 
tions used here is a generalization of the results of Wittmer, 
Beckrich et a/.^r^i^^^ who used a equivalent form of ROL the- 
ory to successfully predict universal corrections to Gaussian 
chain statistics for homopolymers in a dense melt. In a diblock 
copolymer melt, the predicted correction to the random-walk 
model for rj^j (q) is given by a function of the form 



cN 



6Q,,{qRg,xoNjA,bB/bA), (24) 



and is smaller than the ideal gas contribution ilij (q) by a pref- 
actorof Here, fi^, is a dimensionless function that, in 

the case of interest here, can only be evaluated numerically, by 
a procedure outlined in App. |A] The corresponding correction 
to the random walk model for F{q) is given by 

<5F(q) = -ciV^e,f],r.i(q)5r!,,(q)iife/(q)ei. (25) 



This is the expression that we use to evaluate the intramolec- 
ular one-loop correction to 5^'^(q). 



B. Direct Correlation Function 

The one-loop theory yields an expression for the apparent 
X-parameter as a sum 



Xa(q) = xo + Axa(q), 



(26) 



in which 



AXa(q) - £.£,12-,l(q)I]fe;(q)f2yl(q), (27) 



E,,(q) ^ - y^fi(^,(q,k_,-k+)Gfe,(k+) 
S(3) 



xl7^^^;(-q,-k_,k+)G,„„(k_] 



(28) 



Here, k+ ee k ± 



and is a three point intramolecular 



correlation function for monomers within blocks i, j and k 
of a Gaussian diblock. The physical content of eq. (l28b for 
is shown schematically by the right diagram of Fig. [T] in 
which two pairs of segments along two chains interact via the 
screened interaction G. 



The integral in eq. 



for Sij(q) is also UV divergent. 



and so yields a UV divergent expressions for Axa(q)- It was 
shown in ref.', however, that the divergence of this quantity 
could be interpreted as a renormalization of the interaction 
parameter It was shown that, in a model with equal A and B 
statistical segment lengths, the one-loop prediction for Xa(q) 
may be written as a sum of the form 



Xa(q) = Xe(A) + ^Xa(q), 



(29) 



in which Xe(A) is a cutoff-dependent effective interaction pa- 
rameter that is independent of both q and N, and i5xa(q) is a 
cutoff-independent correction. The correction (5xa(q), which 
depends upon both q and N, represents a correction to the 
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phenomenology predicted by the RPA, rather than merely a 
correction to the value of the interaction parameter It was 
found that this quantity is given by a function of the form 

SXa{q) = j;^^Y/2SXa{qRg,XN, fA,bB/bA), (30) 

where Sxa is a dimensionless function that we calculate by 
numerical integration. 



C. Self-Consistent Approximation 

The FH theory for S'(q) in the disordered phase is based 
on a self-consistent one-loop, or Hartree, treatment of fluctu- 
ations. Such an approximation is obtained by starting from 
a perturbative one-loop approximation for ^^(q), in which 
(55 (q) is initially expressed as a Fourier integral involving the 
RPA or "bare" correlation function 5*0 (q), and then replacing 
Sq (q) by the self-consistently determined correlation function 
5(q) throughout the expression for 6S{q). The use of this 
approximation is justified for symmetric diblock copolymers 
near the ODT by an analysis given by Brazovskii, who con- 
cluded that the Hartree approximation correctly captures the 
dominant contributions to ^^(q) near the spinodal of any ho- 
mogeneous systems for which Landau theory predicts a sec- 
ond order transition to a periodic state. 

To define an analogous self-consistent approximation for 
the auxiliary field theory, we simply evaluate all of our 
expressions for the one-loop correction dftij{qR, xN) and 
6F{qR, xN) by replacing x by an "apparent" x-parameter, 
Xa that is defined by fitting the actual peak-intensity to the 
RPA, by setting 

cNS-\q*)=2[{xN)s-XaN], (31) 

where {x^)s denotes the RPA spinodal value. We thus eval- 
uate the one-loop corrections using a value of x that is chosen 
so that the underlying approximation for 5(q) has the correct, 
self-consistently determined peak value. In this approxima- 
tion, we have 

5-i(q) = So\ci)+SS-\(iR,XaN). (32) 

Here, the function 6S^^{qR, XaN) represents the one-loop 
correction obtained from the perturbative one-loop theory, af- 
ter subtracting UV divergent terms. In the same approxima- 
tion, the microscopic parameter Xe (which is assumed to de- 
pend on temperature and details of local liquid structures, but 
not on N) is related to Xa by a relationship 

cN 

XaN = XeN- —6S-\q*R,XaN). (33) 

The self-consistent theory can be evaluated numerically by 
calculating the one-loop contribution for a specified sequence 
of values of Xa^, and then using eq. (|33l l to infer correspond- 
ing values of Xe-^- 




(a) 



FIG. 2. Normalized correlation function cNS{q) for symmetric di- 
block copolymer with Xe = and N — 500, as predicted by the 
RPA (dashed line) and ROL (solid) theories. 

IV. RESULTS 

We present the predictions of the ROL theory in this sec- 
tion, and compare them to the RAP and FH theories when- 
ever possible. Since we are interested in the disordered phase 
throughout this work, from now on, we use wavenumbers 
q = |q| instead of vectors q as arguments. 

A. "Ideal" Symmetric Copolymers (x = 0) 

We begin by considering the special case of a model with 
Xo = 0, in which the underlying pair interaction Uij{r) is 
independent of the monomer type indices i and j. This de- 
scribes a homopolymer melt in which we arbitrarily label the 
first fj^N monomers of each chain as A monomers, and the re- 
mainder as B, but in which A and B monomers are physically 
indistinguishable. This situation is approximately realized in 
neutron scattering experiments in which the two blocks of a 
diblock molecule contain protonated and deuterated versions 
of the same monomer. We refer to this as an "ideal" copoly- 
mer by analogy to the description of a mixture of two almost 
indistinguishable molecular species as an "ideal" solution. 

We further focus our attention on the case of symmetric 
diblock copolymers, with = /s = 1/2 and bA ~ 6b, and 
in which the two blocks are also structurally identical. We 
show in App. |B] that, in this case, Xa (q) — for all q when 
Xo = 0. The inverse correlation function of an incompressible 
melt of such polymers thus is given exactly by 

cNS-\q) = F{q), (34) 

where F{q) is the combination of intramolecular correlation 
functions defined in eq. ( fT4l i. In this case, the deviation of 
S^^{q) from Leibler's RPA prediction for ideal diblocks is 
thus entirely a result of a deviation from Gaussian single-chain 
statistics. 

Fig. |2(a)| shows predictions for the normalized structure 
factor S{q)/cN for a melt of symmetric diblock copolymers 
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FIG. 3. Normalized one-loop correction N^^^SF{q) to the inverse 
correlation function for symmetric diblock copolymers with Xe = 0. 



with N — 500 and Xe — 0, as predicted by original RPA 
theory (dashed line) and by the ROL theory (solid line). The 
one-loop corrections to the RPA cause a slight increase in peak 
position g*, to a value slightly above the RPA value qq, and 
cause a slight enhancement in peak intensity S{q*). Both of 
these trends are opposite to those found near the ODT, where 
composition fluctuations tend to decrease q* and suppress the 
peak intensity. For large q values, qRg > 5, the correction to 
the RPA becomes negative, and is found to decrease asymptot- 
ically as q~^ for qRg ^ 1. The behavior in this high-q regime 
is identical to that found by the Strasbourg group^^-^^ for the 
high-(7 corrections to the intramolecular correlation function 
in a homopolymer melt. 

The predicted corrections to the RPA are proportional to 
7V~^/^, and so would be smaller than shown in Fig. |2(a)| l for 
longer chains. In Fig. IIV Al we show the one-loop correc- 
tion SF{q) to the random walk model for F{q), multiplied 
by N^/'^ to obtain a universal function of qRg that is inde- 
pendent of N. Three different wavenumber regimes are visi- 
ble in this representation: For small wavenumbers, qRg < 1, 
both Folq) and 6F{q) are positive, and both diverge as l/q^ 
as q ^ 0. We show in Sec. IV that F{q) cx l/{qRAB)^ 
for qRg <C 1, where R\g is the mean-squared separation 
between the centers of mass of the A and B blocks. A pos- 
itive correction to F{q) in this limit thus corresponds to a 
negative correction to the random walk prediction for Rab- 
At very high wavenumbers, qRg ^ 5, we find SF{q) ^ 
N-^''^{qRgf and Fo{q) - {q^Rgf, which leads to the be- 
havior cNdS{q) ~N^^/'^{qRg)^^ noted above. In the 
intermediate range 1 < qRg < 5 surrounding qoRg ~ 2, 
however 6F{q) is negative and decreases with increasing q, 
leading to an enhancement in S{q) and a slight positive shift 
in q*. 

The ROL theory predictions of a value of q* slightly higher 
than the RPA prediction, and of a value of Rab slightly less 
than the random walk prediction, are meaningful only if ac- 
companied by a clear definition for the value of the statistical 
segment length b used in the random-walk model to which 
the ROL theory is compared. The way the ROL theory is de- 
rived requires that the statistical segment length in the RPA 



theory be defined to be that of a hypothetical system of in- 
finite chains, by defining to be the N ^ oo limit of the 
ratio 6Rg/N in a homologous series of melts of increasing 
chain length N. This is the definition used by the Strasbourg 
group to compare their ROL predictions for the intramolecu- 
lar correlation function in homopolymer melts to the results 
of extensive computer simulations i-^^'-^^-"*^ 



B. Peak Intensity (Symmetric Copolymers, x / 0) 

We now consider how the value of the maximum S{q*) 
in the structure factor for a symmetric diblock copolymer 
evolves with increasing x- Fig- |4]compares predictions of the 
RPA, FH, and ROL theories for the normalized inverse inten- 
sity cNS^^{q'*) vs. XeN for symmetric diblock copolymers 
with N = 1000 and 10000. This corresponds roughly to the 
range explored by most experiments. The straight dashed line 
is the RPA prediction, in which S^^{q*) is a linear function 
that vanishes (or S{q*) diverges) at XeN ~ 10.495. Both the 
FH and ROL theories predict deviations from the RPA that 
decrease as N^^/"^ with increasing chain length, but that re- 
main significant at experimentally relevant chain lengths. The 
dots are the FH theory predictions for the values of x^ at the 
fluctuation-induced first-order transition. The ROL auxiliary 
field theory does not yet provide a prediction for the this tran- 
sition. 

The FH theory and the self-consistent ROL theory predict 
similar results at large values of Xe-^, where the deviations 
from the RPA are dominated by the effects of strong com- 
position fluctuations with wavenumbers (7 ~ This is the 
regime that the FH theory was designed to describe, and in 
which it has been most successful as a description of experi- 
mentally measured scattering intensities in nearly symmetric 
diblock copolymers. ^'Siiirilii^^ Note that deviations from the 
RPA predicted by both FH and ROL theories remain signifi- 
cant even far from the predicted ODT, particularly for the FH 
theory at iV = 1000, that the differences between the FH and 
ROL predictions increase as XeN decreases. The FH theory 
always predicts a positive correction to the RPA expression 
for S''^{q*), or a decrease of S{q*). The ROL theory predicts 
a positive correction to S'~^(g) (a suppression of scattering) 
for XeN > 6 but a negative correction (an enhancement of 
S{q*)) at smaller values of XeN. This allows the ROL pre- 
dictions to interpolate smoothly between the behavior found 
at X = 0, where the theory yields an enhancement of S{q*), 
and behavior similar to that of the FH theory near the ODT. 

Fig. |5] shows a more careful comparison of the asymptotic 
behavior of the ROL in the limit of strong scattering. There, 
we show a log-log plot of the predicted one-loop correction 
N^/^cNSS-'^{qo) as a function of l/[{xN)s - XaN]. The 
square-root divergence predicted by the FH theory, given by 
eq. (|5]l, is shown by the straight line. The convergence of 
ROL towards this line confirms the conclusion of ref.^^ that 
the FH and the self-consistent ROL exhibit the same asymp- 
totic behavior near the spinodal. The differences between the 
two theories become significant for XaN < 9.5. 

The first, and most influential, quantitative comparison be- 
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FIG. 4. Self-consistently calculated inverse peak intensity versus 
XeA'^ for symmetric diblock copolymer with iV = 1000 (a) and 
10000 (b). Predictions of RPA, the FH theory, and the ROL theo- 
ries are shown. The dots marks the ODTs predicted using the FH 
theory. 
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FIG. 5. Correction to (go) predicted by the FH theory (dashed) 
and by the ROL theory (solid) on log-log scale. The two theories 
exhibit the same dominant asymptotic behavior as Xa approaches 

Xs- 



predictions of the FH and ROL theories in the relevant range 
of values of xN, however, using either method of estimat- 
ing N, we find that these theories yield very similar results 
for such high values of N over the limited range of values 
of probed in these experiments. We thus doubt that it is 
possible to distinguish between the two theories by compar- 
ing only their predictions for the disordered state scattering 
intensities repotted in this study. Quantitative comparisons to 
lower molecular weight samples over a wider range of values 
of xN would highlight the differences, as would comparison 
to computer simulations. 



C. Peak Wavenumber (Symmetric Copolymers, x / 0) 



tween the FH theory and experiments was given by Bates 
et al.,^'^ for melts of relatively high molecular weight nearly 
symmetric (/a — 0.55) poly(ethyl-propylene-b-ethyl ethy- 
lene) (PEP-PEE) diblock copolymer. These authors con- 
cluded that the FH theory gave a quantitatively accurate de- 
scription of both the temperature dependence of the peak in- 
tensity near the ODT for systems of fixed molecular weights, 
and of the relationship between the ODT temperatures of sev- 
eral samples of different molecular weights, using a single 
function x(r) = A/T + B. In the original study, N was 
estimated by fitting q*{T) to the RPA prediction, thus assum- 
ing that q*Rg was constant. This gave a value of N that varied 
over the range 6000 < N < 11000 for the most heavily stud- 
ied sample over the experimental temperature range, which 
was found to correspond to a range 10.6 < x^ < 12.5. Us- 
ing data reported in a later work, ' ' we have also estimated N 
by using independently determined statistical segment lengths 
for homopolymers. This procedure yields noticeably lower 
values of N: 3400 < N < 4300. This difference suggests 
that the near perfect agreement between the FH theory and the 
experimental results for both the scattering in the disordered 
phase and the ODT may be partly fortuitous. Upon comparing 



The ROL theory predicts a peak wavenumber q* for sym- 
metric diblock copolymers that decreases slightly with in- 
creasing Xe^, with a more rapid decrease near the transi- 
tion. Fig. |6] shows predictions of the normalized peak posi- 
tion q* /qo vs. x^ for several different chain lengths, where 
qa = 1.946/i?g is the RPA prediction.^ 

Sub-figure (a) shows the peak wavenumber vs. XaN. Since 
the definition of Xa in eq. dSTI ) is directly related to the peak 
intensity, this is essentially a plot of peak wavenumber versus 
a measure of peak intensity. Sub-figure (b) instead shows the 
peak position vs. Xe^- Because Xe is expected to be a func- 
tion of temperature alone, and nearly linear in 1 /T, sub-figure 
(b) shows the predicted behavior of an experimental plot of q* 
vs. N for a sequence of polymers at the same temperature, 
or of q* vs. 1/T in a system for which the statistical segment 
lengths have negligible intrinsic temperature dependence (or 
for which the data has been plotted in a way that corrects for 
this effect). 

The predicted value of q* is slightly larger than go 31 
XeN — 0, as already noted in Sec. IIV Al and less than go ne^" 
the ODT. The dots indicate the FH prediction of the value of 
XeN at the first-order ODT. (The ROL theory cannot be used 
for this purpose, because it has not been used to predict the 
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FIG. 6. Shift in tiie peak position of correlation function versus (a) 
XaN (b) and Xe-^ for symmetric diblocks witii A'^ = 100, 1000 and 
10000, respectively. The sub-figure (b) is prepared using the self- 
consistent approximation discussed in Sec. IIIICI The dots in (b) are 
the FH theory predictions for (xN)odt- 



ODT). The ROL theory predicts a decrease of 10% - 20% in 
q* at the FH ODT, relative to the values at — 0^ over a 
surprisingly wide range of molecular weights. This is roughly 
consistent with the range of values found in simulations on 
values of TV ^ 100 or less and in experiment on systems with 
much higher values of iV: iV ~ 10^ - io4.5.6.8,iQ-i3,i5,i6,i9-2i. 

Quantitative comparisons to both simulation and experimental 
data will be presented elsewhere. 

The observation that q* tends to decrease with increasing 
xN near the ODT has sometimes been described, somewhat 
loosely, as a result of "chain stretching" induced by strong 
fluctuations near the ODT.'^ This description implies that the 
shift in q* is a result of change in intramolecular correla- 
tions. Because the ROL theory provides predictions of both 
intramolecular and collective coiTelation functions, it is pos- 
sible for us to ask whether this is the correct interpretation. 
Our discussion of this question is based upon the generalized 
Ornstein-Zemicke equation, eq. (fTTT i. in which cNS^^{q) 
is given as the sum of a term F{q) that depends only upon 
single-chain correlation functions and a term —2xa{q)N that 
is sensitive to inter-molecular coiTelations. If the shift of q* 
were primarily a result of changes in the intramolecular cor- 
relations, we would expect it to be well approximated by a 
modified RPA theory in which Leibler's approximation for 
F{q), which assumes Gaussian chain statistics, is replaced by 
the function that is obtained from the predicted, slightly non- 
Gaussian intra-molecular correlation functions, but in which 
we still approximate Xa{q) by a wavenumber independent 
value. 

To clarify the origin of the dependence of q* upon xiV, Fig. 
|7] shows the g-dependence of the predicted inverse structure 
factor S~^{q) and its components, for a system with N = 500 
and XeN = 8.0. This figure shows the ROL prediction for 
S~^{q), the RPA prediction SQ^{q), and the one-loop con- 
tributions SF{q) and -2Sxa{q)N to F{q) and -2xa{q)N, 
respectively. For q near qo, SF{q) decreases with increasing 
q, while —2xa{q)N increases. Ignoring the q dependence of 
~Xa{q), as suggested above, would thus yield a value of q* 




XeN = 8.0 

12 3 4 

qRg 

FIG. 7. Wavenumber dependence of inter- and intra- molecular com- 
ponents to the total correlation functions from the one loop predic- 
tion, evaluated for N — 500 and at x>i^ ~ 8.0. Following eq. il3\ . 
to display contributions from each component, the inverse of S{q) 
is shown. cNSq^: Leibler's mean field theory prediction. F{q): 
the intramolecular part with one-loop correction. ~2xa{q)N: the 
intermolecular part with one-loop correction. 



slightly larger than qo. This is indeed what is found when 
xN = 0, where Sxa{q) — 0. The sign of this effect is in- 
dependent of xN: The contribution to 6F{q) that arises from 
deviations from Gaussian chain statistics is always a decreas- 
ing function of q that would, by itself, always tend to increase 
q*. The fact that q* drops to values lower than go ns^r the 
ODT is a direct result of the fact that -~2Nxa{q) increases 
with increasing q, and that this effect dominates near the ODT. 
We conclude that the decrease in q* with increasing xiV is a 
result of intermolecular, rather than intramolecular, correla- 
tions. Similar conclusions were reached by Binder and Fried 
on the basis of Monte Carlo simulations.^^ 

Fig. [8]shows values of the ROL predictions to SF{q*) and 
-~2xa{q*)N at the peak wavenumber q* as functions of XaN , 
for a melt of symmetric copolymers with N — 1000. The 
intramolecular contribution 6F{q*) is always negative, and 
thus enhances the peak intensity S{q*), while —Sxa{q*)N 
is positive, and thus suppresses the peak intensity. The pre- 
dicted crossover from negative values of SS^^{q*) (enhanced 
scattering) at low values of xN to positive values (suppressed 
scattering) closer to the ODT is caused by a change in the 
relative importance of these two contributions with increasing 
xN. 



D. Peak Intensity (Asymmetric Copolymers) 

The ROL theory may also be used to describe asymmetiic 
copolymers. We have shown elsewhere^^ that both the full 
MFEH theory and the ROL theory yield one-loop contribu- 
tions to S^^{q) that, in general exhibit behavior near a spin- 
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FIG. 8. The corrections to the inverse of peak intensity evalu- 
ated using the ROL theory, alongside with the component contribu- 
tions from inter and intra molecular correlations, respectively, versus 
Xa^ ■, for TV = 1000. The two components have opposite signs near 
the mean field spinodal. 




odal that can be described by an expansion of the form 



1 



iVl/2 



B 



(35) 



where the coefficients A, B and C depend on copolymer com- 
position Ja and on the ratio of statistical segment lengths. 

The coefficient A of the highest order 1/r divergence can 
be shown to vanish in the special case of a completely sym- 
metric copolymer (/yi = /b and 6^ = bs) considered by FH, 
as a result of the symmetry between the two blocks. These au- 
thors took advantage of this symmetry by dropping a contribu- 
tion to the MFEH one-loop theory (a Fourier integral) that is 
zero in this special case, but that actually yields the dominant 
l/r divergence for asymmetric copolymers. It is straightfor- 
ward to show that the coefficients A and B for asymmetric 
copolymers are opposite in sign: The coefficient B that is re- 
tained in the FH theory is positive, and thus tends to increase 
S~^{q) or decrease S{q), while A is negative for /a 7^ 0, and 
thus tends to increase S{q) near the spinodal. 

Fig. (|9]l shows results of the ROL theory for the normalized 
one-loop contribution N^/^cNSS~^{qo) for modestly asym- 
metric and symmetric diblock copolymers with 6a — bs, for 
systems with /a = 0.35, 0.40, 0.45, and 0.5. In the case 
of a symmetric polymer, the quantity N^/^cNSS~^{qo) in- 
creases monotonically with increasing XaN, as also shown 
Fig. |5] For the most asymmetric case of /a = 0.35, the 
one-loop contribution is clearly a non-monotonic function of 
XaN, which increases with increasing XaN far from the spin- 
odal, but then decreases rapidly near the spinodal. This is a 
result of the competition between the A/t term with ^ < 0, 
and the S/Vr term, with _B > 0, in eq. ( |35] |. 

For asymmetric copolymers, SCFT predicts a first-order 
ODT. The two vertical dotted lines in each figure are SCFT 
predictions for the values of xiV at the spinodal (the larger 
value) and at the ODT (the smaller value). For /a = 0.5, these 
values are equal. Note that, for the most asymmetric copoly- 



FIG. 9. SS-^{qo) versus Xa for diblock copolymer with &a = &s, 
at various block fractions: / = 0.35, 0.40, 0.45 and 0.50. Solid 
lines are fitted results using eq. ( IC2) . The dash-dotted and dotted 
perpendicular lines are ODTs calculated from SCFT and the mean 
field spinodal, respectively. 
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FIG. 10. Self-consistently calculated inverse peak intensity versus 
XeN for diblock copolymer with / — 0.35 and N = 1000. Predic- 
tions of RPA, and the ROL theories are shown as dashed and solid 
lines, respectively. ODT calculated from the SCFT is also marked. 



mer, with /a — 0.35, the maximum in N^/^cNSS^^{qQ) ap- 
pears at a value of xN significantly below the SCFT ODT 
value. The theory thus suggests that strong fluctuation effects 
that are qualitatively different from those found for symmet- 
ric copolymers could appear within the disordered phase of 
asymmetric copolymers. 

In the results shown in Fig. (|9]l, the input parameter xN 
may be calculated using either the SCFT interaction param- 
eter Xe, to obtain a peiturbative form of the ROL theory, or 
using the apparent interaction parameter Xa, to obtain a self- 
consistent Hartree approximation. The distinction between 
these two approaches is unimportant far from the ODT, where 
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both yield small corrections to the RPA, but become important 
near the ODT, where the predicted corrections become large. 
Brazovskii argued that the Hartree approximation should be 
accurate near the crystallization transition in a class of systems 
with a scalar order parameter field that exhibits strong fluctu- 
ations at a nonzero wavenumber q*, and in which a Taylor 
expansion of the effective Hamiltonian has a vanishing cubic 
term (leading to a second-order transition in Landau theory) 
and a small quartic term. His analysis justifies the use of a 
Hartree approximation for symmetric copolymers, for which 
the cubic term in the Landau theory must vanish by symmetry, 
in the limit N oo, in which the effects of the quartic anhar- 
monicity become small. It does not, however, justify the use 
of a Hartree approximation for strongly asymmetric copoly- 
mers. 

Fig. [TO]shows the results obtained if (despite the above dis- 
cussion) we use the self-consistent (Hartree) ROL theory to 
calculate S{q*) for asymmetric diblocks with / — 0.35 and 
N = 1000. For sufficiently small values of Xe^, the behav- 
ior is similar to that shown in Fig. |4]for symmetric diblocks. 
For XeN close to 12, however, we see an unphysical turn- 
ing point, leading to a predicted value of cNS^^{q*) that is a 
multivalued function of Xe^- The behavior is similar to that 
obtained in previous attempts to use a Hartree approximation 
to describe the peak scattering intensity at g = in off-critical 
polymer blends. ^"^'^^ In the example shown here, the turning 
point occurs at a value of Xe^ that is below both the SCFT 
prediction of the ODT for this value of Ja and below the FH 
prediction of {xeN)c — 14.6 for the ODT of a symmetric 
copolymers with N — 1000. Both of these values are pre- 
sumably less than the true value of at the ODT for this 
system. This indicates that the self-consistent ROL is almost 
certainly not an adequate description of moderately asymmet- 
ric copolymers near the ODT. We do, however, expect either 
form of ROL to accurately describe small corrections to the 
RPA fai-ther from the ODT. 



E. Molecule and Block Dimensions 
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FIG. 11.x dependence of the fractional change in the radius of gy- 
ration of one block (-Rg,A) ™d of the whole molecule and the 
center of mass distance of two blocks {B?ab)' for ^ symmetric di- 
block. The curve labelled with represents the corresponding 

fractional change iV^'^^((go/q*)^ — 1) in the length scale (g*)~^. 
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FIG. 12. The change of mean squared distance of center of mass 
of two blocks (Rab^ Isft) and the radius of gyration of one block 
(Rg^A^ right), for block fractions / = 0.1, 0.2, 0.9. The abscissa 
is xN normalized by the mean field spinodal. 



Previous theoretical and simulation studies have character- 
ized the statistics of individual chains primarily by examining 
how the radii of gyration of entire chains and of the A and 
B blocks change with x^- Let Rg be the radius of gyration 
of an entire copolymer, and let Rg^i be the radius of gyration 
of the i block (i —A or B). Let Rab '^^e mean-squared 
distance between the centers of mass of the A and B blocks. 
These quantities are related to each other by 



where we have 



r2 



L 



R^ 



AB 



and an analogous expression for L^g- Using these expres- 
sions, it is straightforward to show that F{q) diverges like 
l/q^ as (7 — > and that 



IaR^ 



fsRl^B + fAfsRAB- 



(36) 



lim F^Vg) 



UaJb^abI ) 



(38) 



They, respectively, are also related to the low-q behavior of 
the intramolecular correlation functions Vlij [q), in the Guinier 
regime qR ^ 1, by relations of the form 



lim ^ij{q) ~ cNfifj 



(37) 



in the same regime. 

Fig. [TT]shows the x dependence of the one-loop corrections 



to Rl, 



i?g A and Rab foi" ^ symmetric diblock copolymer 



Because the predicted fractional changes are all proportional 
to 7V~^/^, we show fractional changes multiplied by N-^^^ to 
obtain quantities that are independent of chain length. 
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At X = 0, the calculated corrections ^, R\g and i?^ are 
all negative. The predicted change in the overall radius of gy- 
ration, SUyR^ = -1.42iV"i/2 at = 0, is the same as 
that found previously for monodisperse homopolymers. ^'^ As 
already noted, these predicted deviations in chain dimensions 
are all defined as deviations from the predictions of a random 
walk model in which the statistical segment length is defined 
to be that of a hypothetical system of infinite chains. Because 
the deviation is negative and decreases as N^^^^ with increas- 
ing N, it implies that longer chains are slightly more swollen 
than shorter ones. The reasons for this A^-dependent deviation 
from random walk statistics has been discussed previously by 
Beckrich et a/.^- 

As X increases, the mean-squared distance R\g between 
the block centers of mass increases, while the block radii of 
gyration decrease slightly. This qualitative trend has been pre- 
viously predicted^' and observed in simulations.^- The phe- 
nomena is consistent with the idea that the two blocks avoid 
increasingly unfavorable AB contacts by shrinking and by 
moving away from one another The deviation of the total 
radius of gyration from the random walk value is negative 
and almost independent of XaN for XaN < 9, but increases 
rapidly with XaN near the spinodal. The predicted corrections 
to all these quantities diverge like t~^/^ close the spinodal, as 
does the correction to the inverse peak intensity. These di- 
vergences are pre-empted by the appearance of a first-order 
transition. For N = 500, the predicted fractional changes in 
Rg, Rab, Rg,A between = and xN = 10 are 2.5%, 
9.4% and -5.5%, respectively. 

For comparison. Fig. [TT]also shows the fractional change 
in the squared peak wavelength 1 / (g*)^ as a function of XaN, 
which increases monotonically with increasing x^. Note that 
the fractional change in l/q^ is predicted to increase notably 
more rapidly with increasing XaN than either 
agreement with simulation results. This comparison is enough 
to show that it is impossible to quantitatively account for the 
shift in value of the nonzero peak wavenumber q* by examin- 
ing overall chain dimensions or (equivalently) the low-g be- 
havior of the single-chain correlation function. The above 
analysis of the behavior of F{q) for q near q* indicates, more- 
over, that any attempt to interpret the shift in q* as a direct re- 
sult of change in intramolecular correlations is actually quali- 
tatively wrong. 

Fig. [12] shows how the same measures of chain and block 
sizes depend on copolymer composition /a in melts of asym- 
metric copolymers. Fig. |12(a)| shows the results for 5R\g 
and Fig. |12(b)| for 5R^g ^. In both plots, the x axis is the ratio 
of XaN to its spinodal value. 5R\g shows a relatively weak 
dependence on /. For highly asymmetric copolymers, the ra- 
dius of gyration of the minority block shrinks much more than 
the majority block, and the majority block increases with in- 
creasing XaN near the spinodal. 

V. CONCLUSIONS 

We have presented quantitative predictions of the ROL the- 
ory for corrections to the RPA description of correlations in 



disordered diblock copolymer melts. The theory is expected 
to have a wider range of validity than the earlier FH theory 
and its descendants, and makes independent predictions about 
intramolecular and collective correlation functions. 

When applied to symmetric diblock copolymers, the ROL 
theory yields predictions that are similar to those of the FH 
and BF theories near the ODT, but that are qualitatively dif- 
ferent in the limit of low x^ values. The ROL theory predicts 
a depression of the peak scattering intensity S{q*) near the 
ODT but a slight enhancement in S{q*) for xN < 6, and de- 
crease in q* near the ODT but a slight increase at small values 
of xiV, relative to RPA predictions. In the limit xN — 0, these 
deviations from RPA predictions are direct results of the de- 
viations from Gaussian single-chain statistics in a dense poly- 
mer liquid that have been studied previously by the Strasbourg 
group. Near the ODT, corrections to the RPA are instead 
dominated by corrections to the apparent interaction parame- 
ter Xa(g)- 

Like the BF theory, the ROL theory predicts a monotonic 
decrease in q* with increasing xiV. The ROL theory also pre- 
dicts a shrinkage of individual blocks, and an increase in the 
distance between the A and B blocks, in agreement with sim- 
ulation results. We have tried to further clarify the physical 
origin of the shift in q* by comparing the results of the full 
theory to those of a simplified theory in which we take into 
account the predicted deviations from single-chain statistics, 
via the function F{q), but approximate the apparent interac- 
tion parameter Xa (?) by a parameter that is independent of q. 
In this simplified theory, the effect of deviations from Gaus- 
sian statistics alone would actually be to increase q* above the 
value predicted by the RPA, as is found in the regime xN < 6 
in which this intramolecular effect dominates the correction to 
S{q). The decrease in q* with increasing x^ is instead found 
to be a result of the (/-dependence of the apparent interaction 
parameter, and is thus best understood as a result of changes 
in intermolecular rather than intramolecular correlations. 

When applied to asymmetric diblock copolymers, the ROL 
theory predicts a correction to the RPA prediction for S^^{q*) 
that increases with increasing xN far from the spinodal transi- 
tion, but that reaches a maximum and then begins to decrease 
very near the spinodal. Applying the "Hartree" approximation 
that underlies the Brazovskiiand FH theories to asymmetric 
polymers yields strikingly unphysical results near the ODT. 

A comparison of the predictions of the ROL and FH theo- 
ries for symmetric copolymers over the range of parameters 
examined in the neutron scattering data of Bates, Rosedale 
and Fredrickson suggests that it would be difficult to dis- 
criminate between these two theories on the basis of these 
experiments, which were carried out over a limited range of 
values of xN in a model system of unusually high molecu- 
lar weight. We will compare both theories to the results of 
computer simulations of shorter polymers over a wide range 
of values of xN in a subsequent paper 
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Appendix A: Renormalization Procedure 

The results presented in this paper are all obtained by ap- 
plying the renormalization scheme that was outlined in ref. ' to 
the Fourier integrals that appear in eqs. (l2ll and (I2FI . In this 
scheme, the UV divergent parts of the Fourier integrals that 
contribute to S^^ (q) are all identified with renormalization of 
the interaction and statistical segment length parameters, and 
thus implicitly absorbed into changes in the values of these pa- 
rameters, while the UV convergent parts are identified as cor- 
rections to the RPA. To extract the UV convergent part of such 
integrals, we have numerically evaluated each integral with re- 
spect to a wavevector k over a domain |k| < A, for AR ^ 1, 
and then subtracted off the analytic result obtained in ref. ' for 
the UV divergent part of the corresponding integral. In each 
case, the resulting difference is found to be approximately in- 
dependent of A for AR 1, and to approach a finite limit 
as AR 00, confirming our analytic calculations of the UV 
divergent contributions. Our final estimate for the UV conver- 
gent part of each such integral is obtained by repeating this 
procedure for several values of A and then numerically ex- 
trapolating to A = 00. The integrals with respect to k in eqs. 
(l2Tl i and ( |28] | were calculated with a 2D Rhomberg integration 
algorithm,'*'* using the wavenumber A; = |k| and the cosine of 
the angle between k and q as coordinates. 

A similar procedure was used in ref^^ to calculate the one- 
loop corrections to the free energy density and to the q ^ 
limit of Xa (q) in a polymer blend. These calculations only 
required one-dimensional integrals with respect to |k|. 

Intramolecular Correlations It was shown in ref. ' that the 
result of eq. dHJ for j (q) can be expressed as a sum 



17,, (q) = r!,,(q) + Wff^ + <5r!,,(q), 



(Al) 



in which is a UV divergent contribution that was found 
to be of the form 



dbk 



Sbk{A), 



where 



Sbk{A) = llA/inH) 



(A2) 



(A3) 



is a UV divergent shift in the value of the statistical segment 
length bk- Here, / — JaIa + JbIb, where li — v/bf. The 
function 17^ (q; b) denotes the random-walk model for fiy (q) 
for a chain with specified statistical segment lengths bA and 
bs- The UV divergent term was thus interpreted as the first 
term in an expansion of the intramolecular correlation func- 
tion rtij (q, b + 6b) for a random walk with renormalized sta- 
tistical segment lengths. The UV convergent term (5f2y (q) is 
calculated by subtracting the result of eqs. (IA2l i- (IA3l l for the 
UV-divergent part of the integral from a numerical result for 
the full integral. 

Direct Correlations The one-loop correction to Xa (q) is ob- 
tained from eqs. (|27] | and (|28] |. It was shown in ref^^ that this 
quantity can be expressed as a sum of the form 



Xa(q) = Xe(A) 



HjqR) 

N 



A + (5x(q), 



(A4) 



with a cutoff-dependent effective interaction parameter Xe (A) 
of the form 



Xe(A) =xo + ^A3 + BxoA. 



(A5) 



The coefficients A and B are independent of A, N and q, 
but depend upon the parameters bA, ^s, v and /a- The co- 
efficient H{qR) depends upon the normalized wavenumber 
qi? as well as bA, bs, v and Ja- Explicit expressions for 
A, B, and H are given in ref.*. We argued there that the 
A-dependent terms in eq. ( IA5l l were naturally interpreted as 
trivial corrections to the value of the local interaction param- 
eter, because they are found to be independent of N and q, 
but do depend on parameters that reflect the local liquid cor- 
relations of this simple model. By the same reasoning, the 
contribution H{qR)A/N cannot be interpreted as a simple 
renormalization of x, because it depends on both q and N. It 
was shown, however, that this term could be attributed to cor- 
rections to the RPA that arise from: i) perturbations in local 
liquid structure near chain ends and near the junction between 
the A and B blocks, and ii) a gradient-squared contribution 
to the excess free energy functional. It was also found, how- 
ever, that the coefficients H{qR) and A both vanish in the 
one-loop approximation in the special case bA = bs- The re- 
sults presented in this paper are all obtained for systems with 
bA ~ bs, for which the only UV divergent contribution to 
Xa(q) is a constant BxqA. The UV-convergent contribution 
Sxa (q) has thus been calculated by subtracting this contribu- 
tion from a numerical result for the full one-loop correction 
AXa(q,A). 



Appendix B: Ideal Symmetric Copolymers 

Here, we show that Xa (q) = for a completely symmet- 
ric diblock copolymer with xo = 0. We first give a general 
symmetry argument as to why this must be true in this special 
case, and then confirm explicitly that this is consistent with 
results obtained from the one-loop theory. 

General Symmetry Argument: Consider the relationship be- 
tween the correlation function matrix S'y (q) in a liquid of 
diblock copolymers containing physically identical A and B 
blocks of equal lengths and the correlation function S{q) in a 
corresponding homopolymer liquid, in which we make no dis- 
tinction between the two blocks. In the homopolymer melt. 



5(q) -r!(q)+i/(q), 



(Bl) 



where ri(q) and H{q) are the intra- and inter-molecular pair 
correlation functions, respectively, in the liquid with no dif- 
ferential labelling. If each molecule in this liquid is instead 
treated as a diblock, we obtain a matrix of correlation func- 
tions 



^»(q) 



(B2) 



where ftij (q) is a matrix of intramolecular correlation func- 
tions, and (q) /4 is an intermolecular contribution, in which 
H{q) is the function defined for a one-component liquid in 
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eq. ( IBll i. The scalar intramolecular correlation function ^{q) 
of eq. (IB II) is given by the sum r2(q) = J^ij ^iji'i)- Th^ 
intermolecular contribution to Sij (q) is simply H{q) /4 in eq. 
(IB2b because each inter-molecular pair in the one-component 
liquid has an equal probability of being labelled AA, BB, AB 
or BA by a process in which one end one end of each chain is 
labelled A and the other B at random. 

Because the matrices 5*^ (q) and ilij (r) must be invariant 
under an exchange of the labels A and B, the eigenvectors of 
Sij (q) and ilij (q) must thus be proportional to the even and 
odd vectors 



5^(1,1), 

£^(1,-1). 



(B3) 



It follows that Sij (q) may be expressed as a sum 

5y (q) = S+iq)S,Sj/2 + 5_(q)e,£,/2, (B4) 

in which S'+(q) and S'-(q) are its eigenvalues. Let i^+(q) 
and fi-(q) be the corresponding eigenvalues of r2ij(q). By 
contracting eq. (fTOl i for 5, j(q) with e to project onto the odd 
subspace, it is straightforward to show that 



(B5) 



By projecting the Ornstein-Zernicke equation for S~j^{q) 
onto the odd subspace, and using the fact that S'r^(q) = 
f^I^(q), we found 



"^CijiqjeiEj = 0. 



(B6) 



Appendix C: Empirical Approximations 

To facilitate comparison with experimental and simulation 
results, we have developed empirical approximations for some 
of our results for polymers with bj^ ~ bs- These are presented 
in this appendix. 



1. Peali Intensity 

The behavior of cNSS^^{qo) for diblock copolymers with 
bA = bB but modestly asymmetric compositions is reasonably 
well approximated for 0.35 < / < 0.65 by a function 



cNSS~\qo) 



(CI) 



1 



"iVl/2 



with composition dependent coefficients 

1 



dif) 



a(/) = 0.652 - 0.799 f 



b{f) 



-18.2 ( --/ 



c(/) = 2.50 + 0.857 (^-/ 



d{f) 



-7.58 - 1460 (--/ 



(C2) 



The expression is symmetric with respect to / 1 — /, and 
the coefficient b{f) of the strongest divergence vanishes at 
/ = 1/2. The quality of the fits is illustrated in Fig. |9] 



It follows immediately from eq. (flSl l that Xa(q) = in the 
incompressible limit. 

One-Loop Theory: To calculate 5xa(q) in the one-loop the- 
ory, we start from eqs. (l27t and (|28) for SCij (q). When xo = 
0, eq. ( fT9] l for Gij (q) reduces to a scalar G(q) — 1 /i^(q). As 
a result, the summation over fc, m and n in eq. (|28] | can be 
carried out separately. By defining 

f],(q,k) ^ ^ar/(q)f](l(q,k_,-k+), (B7) 

kirn 

we may rewrite eq. (l27t as 

5C^M)^\ j f^.(q,k)f],(q,k)G-l(k+)G-l(k_), 

(B8) 

where we have used the fact that f2i(q, k) = f2i(— q, — k). 

For a symmetric diblock copolymer, however, we also have 
f2yi(q, k) = i^s(q, k), as a result of the symmetry under an 
exchange of the labels of the A and B monomers. Using this in 
eq. ( |B8l l yields an expression for SCij (q) that is independent 
of the indices i and j. 



2. Peak Wavenumber 

Let /(x, y, N) denote the quantity cNS^^{q), where x = 
qRg,o and y = xN. The ROL theory prediction for this quan- 
tity is a sum of the form 



f{x,y,N) = fo{x,y) + 



hix, y) 

iVl/2 ■ 



(C3) 



in which /o(x, y) — cNSq ^{qR, x^) is the RPA prediction, 
and /i {x,y)/ N^/"^ is the one loop contribution. 

Let xq be the value of x at which the fo{x,y) is mini- 
mized, so that xq — qoRg.o — 1.946 for a symmetric diblock 
copolymer Because the RPA prediction fo{x, y) is linear in 
y — xN, the minimum with respect to x is independent of y. 

When the minimum of f{x, y, N) with respect to a; = qR 
remains near xq, we may approximate it by Taylor expanding 
both /o {x, y) and /i (x, y) about xq to quadratic order in a; — 
Xq. By minimizing the resulting approximation for / with 
respect to x, we obtain 



Xq 



fiixQ.,y) 



f(;{xo,y) + N-y^fi'(xo,y[ 



■N~^''^. (C4) 
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XaN 

FIG. 13. dependence of the peak position q*, normalized by the 
RPA theory value qo, for N = 100, 1000 and 10000. Solid lines: 
exact results. Dashed lines: estimates using eq. iC4\ . 



For symmetric diblock copolymers, /g'(a;o) = 14.5715. Our 
numerical results for f[{xo;xN) /"(xo;xiV) for sym- 
metric diblock copolymers are accurately fit by the following 
empirical formulas: 



f['{xo,xN) = 



-54.335 + 14.713xiV 

VXsN - xN 
-117.762+ 13.459x7V 



+ 



VXsN-xN 
477.324 - 45. 569x7V 
X.N-xN ■ 



(C5) 



(C6) 



To obtain the self-consistent ROL theory, the parameter x 
should be replaced by the self-consistently determined value 
of XaN throughout the above. The resulting approximation 
for q* is compared to results of the full theory in Fig. [13] 
The quality of the fit evidently improve as N increases, and is 
adequate for most purposes for N > 500. 
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